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Abstract 

The  constant  volume  (CV)  limit  of  pulsed  propulsion  was  explored  theoretically,  where  the  combustion 
chamber  was  approximated  as  being  time-varying  but  spatially  uniform,  while  the  nozzle  flow  was  approximated  as 
being  one  dimensional  but  quasi-steady.  The  CV  limit  was  explored  for  the  isentropic  blow  down  of  a  constant  y 
ideal  gas  for  fixed  expansion  ratios  and  for  variable  expansion  ratios  which  were  adjusted  to  match  the  pressure  ratio 
at  all  times.  The  CV  limit  calculations  were  notable  in  that  all  the  fixed  expansion  ratio  results  could  be  expressed  as 
analytical  solutions.  The  CV  limit  was  then  compared  with  two  relevant  constant  pressure  (CP)  cycles.  Among  the 
several  conclusions  were  that  using  a  fixed  expansion  ratio  nozzle  during  a  CV  blow  down  did  not  exert  more  than  a 
3%  performance  penalty  over  using  a  variable  expansion  ratio  nozzle  as  long  as  the  fixed  expansion  ratio  nozzle  was 
optimized  for  the  blow  down.  Another  major  conclusion  was  that,  while  the  impulse  produced  by  a  CV  device  could 
significantly  exceed  that  of  a  CP  device  operating  at  the  fill  pressure  of  the  CV  device  at  elevated  ambient  pressures 
{e.g.,  1  atm),  the  impulse  produced  by  a  CP  device  could  actually  exceed,  albeit  only  slightly,  the  impulse  produced 
by  a  CV  device  when  the  ambient  pressure  was  near  zero,  such  as  would  occur  in  the  near  vacuum  conditions  in 
space. 

T  -  r/coA*V  ;  dimensionless  time 
Subscripts 

cp  -  pertaining  to  a  constant  pressure  cycle 

e  -  nozzle  exit 

/  -  fill  condition 

Ic  -  pertaining  to  limit  cycle  operation 

0  -  initial  condition  before  start  of  blow  down 

00  -  ambient  condition 

Superscripts 
*  -  throat  condition 

1.  Introduction 

Recent  interest  in  pulsed  detonation  propulsion  has 
spawned  a  number  of  attempts  to  model  the  system 
performance,  the  results  from  which  have  so  far  varied 
widely'.  The  constant  volume  (CV)  limit  of  pulsed 
propulsion,  as  used  here,  refers  to  a  limiting  case  for 
pulsed  propulsion  cycles  which  is  approached  when 
blow  down  times  become  much  longer  than 
characteristic  wave  transit  times  in  the  combustion 
chamber.  Compared  with  how  a  practical  pulsed 
propulsion  device  is  likely  to  operate,  the  CV  limit 
probably  underpredicts  the  thrust  (unrealistically  small 
nozzles),  but  the  specific  impulse  is  probably  close  or 
equal  to  an  upper  bound.  Thus  the  CV  limit  is  a  useful 
reference  case.  Some  analytical  solutions  for  the  CV 
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Nomenclature 

A  -  area 

c  -  speed  of  sound 

Cf  -  thrust  coefficient 

Cp  -  specific  heat  at  constant  pressure 

F  -  thrust 

g(Y)  -  eq.  (4) 

I  -  impulse 

m  -  rate  of  mass  flow 

M  -  Mach  number 

V  -  velocity 

P  -  pressure 

r  -  P/Po 

R  -  gas  constant 

1  -  time 

T  -  temperature 

V  -  volume 

Greek 

£  -  expansion  ratio 

0  -  P/P„ 

y  -  ratio  of  specific  heats 

i  -  1 1  PqCqV  ;  dimensionless  impulse 

p  -  density 

£2,  <P-  eq.  (23) 
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limit  are  explored  here  for  the  isentropic  blow  down  of 
a  constant  y  ideal  gas  at  fixed  expansion  ratios.  The 
results  are  compared  with  the  case  where  the  expansion 
ratio  is  continuously  varied  to  match  the  pressure  ratio 
at  all  times  during  the  blow  down. 


2.  Formulation 

2.1  General 

The  impulse  produced  by  the  unsteady  blow  down 
of  the  combustion  gases  in  the  CV  limit  is  the  integral 
of  the  thrust  F  =  mv^  +  -  P„)A^  over  time,  where 

m  is  the  mass  flow  rate,  is  the  ambient  pressure, 
and  ,  P^,  and  are  the  velocity,  pressure,  and  area, 
respectively,  at  the  exit  plane.  Under  the  transformation 
dt  -  {dt  /  dp)dp  =  -{V I  m)dp  and  the  further 

transformation  r  =  plpQ,  where  V  is  the  combustion 
chamber  volume  and  Pq  is  the  initial  chamber  density 
after  combustion  but  before  blow  down  at  time  ?  =  0 , 
the  total  impulse  I  and  blow  down  time  t  may  be 
expressed  as 

7(0  =  =  PoVf'iVe  +  (Pe  -  P^)K  ,  (1) 

t  =^^dt  =  Pf^  ^dr!  m  (2) 


2.2  Fixed  Nozzles 

For  isentropic  flow,  the  instantaneous  exit  velocity 
is  =  ^2Cp(T -T^)  ,  where  T  and  are  the 

instantaneous  chamber  and  exit  temperatures, 
respectively.  Pressures  and  densities  will  be  related  in  a 
simple  fashion  to  temperatures  according  to 

r/To  =(p/Po)’'‘‘ where  y  is  the 
ratio  of  specific  heats  and  where  the  subscript  “0” 
denotes  conditions  in  the  chamber  at  time  /  =  0 .  With 


Cq  =  -y/yRTo  and  =  y7?  /(y  - 1) ,  the  exit  velocity  can 
be  expressed  as 

^  n1/2 


=  ..Ur-i)/2 


^  2  ^ 


LQ/ 


.r-i 


(3) 

,r-ij  ’  ' 

where  r^  =  pj  p  is  the  ratio  of  the  density  at  the  exit 
plane  to  the  density  in  the  chamber.  For  isentropic 
nozzle  flow,  the  exit  density  ratio  is  related  to  the 
expansion  ratio  e  =  I  A* ,  where  A*  is  the  throat 
area,  by  {see  eq.  5.3  in  ref.  2): 

.  _  SiY)  . 


g(Y)  = 


ry-n 

1/2 

1  2  J 

[y  +  lj 

.r-i 


N(l/2)(y+l)/(y-l) 


(4) 


Equation  (4)  has  two  solutions  for  r<,  corresponding  to 
subsonic  and  supersonic  flow  at  the  exit  plane.  These 
solutions  will  depend  only  on  e  and  y  .  For  supersonic 

flow,  the  density  ratio  will  therefore  remain  constant 
during  blow  down  as  long  as  compression  waves  do  not 
enter  the  nozzle,  and  eq.  (3)  will  be  analytically  inte- 
grable  when  inserted  into  eq.  (1).  Compres¬ 
sion/expansion  waves  will  remain  external  to  the  nozzle 
for  all  r  such  that 


r  >  — 


00 


(5) 


where  pQ  =  P^^  I P^  and  <j)^=P^I  are  the  ratios  of  the 
initial  and  exit  pressures,  respectively,  to  the  ambient 
pressure,  and  where  0'  is  the  critical  value  of  0^  below 
which  a  compression  wave  enters  the  nozzle.  The  criti¬ 
cal  value  of  0^  is  given  by  {see  eqs.  2.48a  and  5.2  of 
ref.  2) 

2y 


¥e 


1-1 


H-- 


y-i-1 


where 


mI=- 


y-l 


/  1  V-' 


-1 


(6) 


(7) 


While  eq.  (5)  holds,  the  flow  will  also  be  choked, 
and  the  throat  conditions  will  be  functions  only  of  y . 
For  a  throat  density  and  velocity  given  by  {see  eq.  2.35a 
of  ref.  2)  p*  =  rpo[2/(y-Hl)f'’"'>  and 

V*  =  -^yRT*  =  Co(2/(y  ,  and  with 

Pe  ~  Pq^I  ’  ^he  instantaneous  mass  flow  and  thrust 
become 


m{r)  =  p  Av  =  PqCqA  g(y)| 
F(r)/PoA*g(y)  =  r^y 

y-1 


y-l 


y-l  , 
+ 


(9) 


1 


where  eq.  (4)  has  also  been  used.  Thus  all  terms  in  eq. 
(1)  will  integrable  analytically. 

Define  i  =  IIpqCqV  and  T  =  tCQA*IV  to  be  a  di¬ 
mensionless  impulse  and  a  dimensionless  time,  respec¬ 
tively.  Then  integrating  eqs.  (1)  and  (2) 
gives 
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1]  _ 

\  ^  1 

1  /  ^ 

[  ^  1 

[r-ij 

[y  +  lj 

1  1 

fy-0 

/i  ,.r-i 

2 

\  J 

1/2 


Mr)  + 


y+1  00 


^  2  ^ 
y-1 


r(r)  = 


sir) 


r  n  \ 


y-1 

V  '  / 


1/2 


b(r). 


(10) 

b(r) 

(11) 


where 


b{r)  =  (1/ r)^'>'~^^'^ .  (12) 

The  quantity  t  /  T  will  be  related  to  the  average  thrust 
F  =  1  /t .  Noting  that  Cq  =  yRTg  and  Pq  =  PqRTq  ,  then 
combining  constants  after  dividing  leads  to 

Yilr  =  FIPaA*=Cp,  (13) 

where  Cp  is  the  average  thrust  coefficient  for  the  blow 
down. 


2.3  Variable  Nozzles 

Fixed  expansion  ratio  nozzles  will  be  optimized  for 
at  most  only  a  single  value  of  r  during  the  blow  down; 
losses  will  occur  at  all  other  values.  To  assess  what  the 
performance  would  be  if  these  losses  were  not  present, 
a  variable  nozzle  can  be  envisioned  where  the  expan¬ 
sion  ratio  is  continuously  adjusted  to  match  the  pressure 
ratio  at  all  times.  The  variable  nozzle  limit  produces  the 
maximum  possible  impulse,  and  is  therefore  also  useful 
as  a  reference  case. 

With  the  exit  pressure  always  matched  to  the  ambi¬ 
ent  pressure,  only  the  exit  velocity  term 

=  ^2Cp(T -Tg)  in  eq.  (1)  will  be  of  concern.  For 

fixed  nozzles,  it  was  found  that  the  term  T^IT  =  r/"* 
was  a  constant,  but  here  the  exit  temperature  is  fixed  by 
the  exit  conditions  P^=P^,  =  Po(l/0o)'^’^ .  and 

=Tq(1/(I>q)^'^~''^'^  .  Therefore  T^/T  will  not  be  con¬ 
stant  in  this  case.  The  dimensionless  impulse  in  this 
case  becomes 


Eq.  (14)  cannot  be  integrated  analytically  but  may  eas¬ 
ily  be  integrated  numerically. 

The  blow  down  time  for  the  variable  nozzle  limit 
will  depend  on  the  mass  flow,  and  calculation  of  the 
mass  flow  will  depend  on  whether  the  exit  area  or  the 
throat  area  is  varied  to  match  the  pressure  ratio.  The 
mass  flow  is  best  calculated  using  the  area  which  is 
being  held  constant.  If  the  exit  area  is  assumed  to  be 
varied,  as  might  approximately  be  the  case  with  an 
automatically  compensating  nozzle  such  as  a  plug  noz¬ 
zle,  the  mass  flow  and  blow  down  time  based  on  a  con¬ 


stant  throat  area  are  still  given  by  eqs.  (8)  and  (11),  re¬ 
spectively.  If  the  throat  area  is  varied,  then  the  mass 
flow  and  blow  down  time  for  a  fixed  exit  area 
are 


bt  =  PeA^v,  =  PoCqA, 


sl/2/ 


y-1 


■Aly 


(15) 


T^('')  =  0o 


llY 


y-1 


1/2  - 

I  (16) 


where  =  IcqA^IV  is  defined  this  time  in  terms  of  the 
exit  area.  Eq.  (16)  also  cannot  be  integrated  analytically 
but  may  easily  be  integrated  numerically.  Inasmuch  as 
shocks  and  expansion  waves  will  not  occur  when  the 
expansion  ratio  is  always  matched  to  the  instantaneous 
pressure  ratio,  eqs.  (14)-(16)  above  for  the  variable 
nozzle  case  will  apply  until  the  expansion  ratio  reduces 
to  unity  and  the  flow  becomes  unchoked,  i.e.,  for  all  r 
such  that 


r  > 


1  2  J 

^00, 

(17) 


2.4  Limit  Cycle  Operation 

In  repetitive  operation,  intake  valves  would  be 
opened  at  some  point  during  the  blow  down  in  order  to 
introduce  fresh  propellants  for  the  next  cycle.  Any  re¬ 
maining  combustion  products  not  yet  expelled  would  be 
pushed  out,  or  purged,  by  the  incoming  fresh  propel¬ 
lants.  If  the  purge  of  the  remaining  combustion  prod¬ 
ucts  can  be  approximated  as  a  constant  pressure  proc¬ 
ess,  and  the  value  of  r  is  large  enough  to  prevent  shocks 
from  entering  the  nozzle,  then  the  above  equations  can 
be  used  to  calculate  the  time  to  purge  the  remaining 
mass  and  the  additional  impulse  produced.  Specifically, 
the  time  required  to  purge  the  remaining  mass  is 
tcp{r)  =  pV / m  =  PqVvI m(r) ,  and  the  additional  im¬ 
pulse  produced  is  /pp(r)  =  F(r)fpp(r) ,  where  F(r)  and 

m{r)  are  given  by  eqs.  (8)  and  (9)  at  a  constant  value 
of  r  during  the  purge.  The  additional  dimensionless 
impulse  and  blow  down  times  can  still  be  expressed  in 
the  forms  of  eqs.  (10)  and  (11),  but  with  the  functions 
a{r)  and  b{r)  given  instead  by 

(18) 

where  the  subscript  "cp”  stands  for  the  these  constant 
pressure  additions.  The  total  impulse  for  the  entire  cy¬ 
cle  can  then  be  expressed,  again  in  the  forms  of  eqs. 
(10)  and  (11),  but  this  time  with  the  functions  a(r)  and 
b{r)  given  by 

a,,(r)  =  a(r)  +  a,/r)  =  l  + 

(19) 


■\)ir 
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bicir)  =  Kr)  +  b^pir)  = 


y  +  1 


(y-l)/2 


-1 


where  the  subscript  “Ic"  stands  for  “limit  cycle.”  For 
limit  cycle  operation  with  variable  nozzles,  the  total 
cycle  impulse  is  given  instead  by  the  sum  of  the  con¬ 
stant  pressure  impulse  with  the  impulse  calculated  by 
eq.  (14),  and  the  total  cycle  blow  down  time  is  the  sum 
of  the  constant  pressure  time  plus  the  blow  down  time 
calculated  by  eq.  (11)  or  eq.  (16).  The  expansion  ratio 
always  matches  the  pressure  ratio  in  the  variable  nozzle 
case,  so  the  exit  density  ratio  must  be  set  to 

r,  =(1/0(1/00)''’'- 


The  density  ratio  r  cannot  be  specified  arbitrarily  if 
limit  cycle  operation  is  to  be  achieved.  Limit  cycle  op¬ 
eration  requires  the  blow  down  to  proceed  at  least  until 
the  chamber  pressure  reaches  the  fill  pressure  at  which 
fresh  propellants  are  introduced.  This  gives 

r  <  {Pj  /  /q)’^’"  ,  where  Pj  is  the  fill  pressure.  The  ine¬ 


quality  sign  indicates  that  the  blow  down  can  proceed 
to  a  pressure  lower  than  the  fill  pressure  if  the  fresh 
reactants  undergo  compression  during  the  fill.  Given 
that  Pq~  Pj  for  CV  heat  addition,  the  ratio  P^/Pq 


could  be  on  the  order  of  1000  times  smaller  for  propel¬ 
lants  in  an  initial  liquid  state  compared  to  an  initial 
gaseous  state. 

For  a  given  density  ratio  r,  the  blow  down  of  the 
remaining  combustion  products  in  a  CP  mode  produces 
a  greater  impulse  than  if  the  blow  down  had  been  per¬ 
mitted  to  proceed  to  a  lower  density  ratio.  Therefore, 
given  a  choice,  the  density  ratio  r  should  be  selected  to 
be  the  maximum  value  possible  consistent  with  limit 
cycle  operation  and  available  injection  pressure,  namely 

r  =  (P//Po)''^ 


2.5  Comparisons  with  Constant  Pressure  Devices 

Comparisons  of  the  dimensionless  impulse  or  the 
thrust  coefficient  with  corresponding  quantities  for  con¬ 
stant  pressure  (CP)  processes  can  be  performed  if  the 
normalization  constants  are  the  same  in  both  cases. 
Here  the  normalization  constants  are  chosen  to  be  those 
of  the  CV  process,  leading  to  the  correction  factors  Q 
and  0  in  the  resulting  expressions  below.  Assuming 
the  expansion  ratio  for  the  CP  process  is  chosen  to 
match  the  pressure  ratio,  the  resulting  dimensionless 
impulse,  blow  down  time,  and  average  thrust  coeffi¬ 
cient  become 


Q 


xl/2 


y-1 


^cp 


0  1 


■^2  g(r) 


Vl-(^/0o)^’'''^'’'  (20) 

(21) 


( 


1/2 


Cf,cp  =  ,  (22) 

''  y-1 0 

where 

Q^CqIc^p  Pol  Pep  (23) 

and  where  =  7,^  /  PqCqV  ,  T,p  =  t^pC^A*  /  V  , 

Pp.cp  =  P I  Po^*  >  the  values  of  Pg ,  Cg  ,  Pg  >  ^'^*3  ^ 
are  those  of  the  constant  volume  process.  The  quantity 
T^p  in  this  context  can  be  interpreted  to  be  the  dimen¬ 
sionless  time  required  by  the  constant  pressure  process 
to  expel  the  same  mass  as  the  initial  mass  of  the  con¬ 
stant  volume  process. 

3.  Results  and  Discussion 

3.1  Effect  of  Thermochemistry 

The  dimensionless  impulse  and  blow  down  time, 
eqs.  (10)-(12),  depend  explicitly  on  r  and  implicitly 
only  on  y ,  (or  e ),  and  0o .  Of  these,  only  y  de¬ 
pends  on  the  thermochemistry,  but  this  dependency  is 
weak.  The  major  influence  of  thermochemistry  comes 
through  the  initial  speed  of  sound  Cg  used  to  normalize 
the  impulse.  Thus  the  specific  impulse  7,  =  iCg  is 
maximized  when  Cg  is  maximized.  This  in  turn  implies 
that  the  optimum  thermochemistry  is  that  which  maxi¬ 
mizes  the  initial  combustion  temperature  Tg  and  mini¬ 
mizes  the  molecular  weight.  The  same  general  guidance 
is  known  to  also  apply  to  constant  pressure  devices. 
Thus  constant  volume  devices  should  optimize  at  ap¬ 
proximately  the  same  mixture  ratios  as  constant  pres¬ 
sure  devices.  Also,  like  constant  pressure  devices,  the 
thrust  is  maximized  by  maximizing  the  initial  chamber 
pressure  and  the  throat  area,  as  implied  by  eq.  (9). 

3.2  Fixed  Nozzles 

Much  can  be  understood  about  the  blow  down  of 
fixed  nozzles  by  examining  conventional  steady  state 
thrust  coefficient  curves  such  as  can  be  found  in  stan¬ 
dard  texts.^  The  thrust  coefficient  is  defined  in  eq.  (13), 
where  Pq  in  the  steady  state  case  is  interpreted  to  be  the 
steady  chamber  pressure  of  a  constant  pressure  device. 
A  set  of  these  curves  for  the  steady  state  case  is  given  in 
Fig.  1  for  y  =  1 .2 .  These  curves  were  generated  by  set¬ 
ting  r  =  l  in  eq.  (18),  substituting  into  eqs.  (10)  and 
(11),  and  then  calculating  Cp  per  eq.  (13).  The  curves 
reproduce  those  given  in  ref.  3.  Dimensionless  isobars 
(curves  of  constant  0g )  initially  increase  with  e  ,  reach 
a  maximum,  then  decrease  to  a  minimum  where  a  shock 
enters  the  nozzle  and  the  above  formulation  is  no  longer 
valid.  The  curve  for  0g  =  °°  reaches  a  maximum  only 
at  e  =  oo  .  A  curve  connecting  the  maxima  indicates  the 
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expansion  ratio  producing  the  maximum  thrust  for  a 
given  isobar. 

The  blow  down  of  a  fixed  nozzle  CV  device  pro¬ 
ceeds  along  a  vertical  path  of  constant  £  between  two 
isobars.  The  thrust  produced  will  be  some  average  be¬ 
tween  the  two  isobars.  Picking  any  two  isobars  in  Fig. 
1,  say  between  0o  =  50  and  0o  =  20,  and  following  the 
vertical  line  between  them  for  various  e  ,  it  can  be  en¬ 
visioned  that  the  average  thrust  will  reach  a  maximum 
neither  at  large  £  nor  at  £  =  1 ,  but  at  some  optimum  £ 
in  between.  However,  the  average  thrust  cannot  be  cal¬ 
culated  directly  from  Fig.  1,  because  the  thrust  coeffi¬ 
cient  is  proportional  to  the  thrust  divided  by  the  cham¬ 
ber  pressure,  not  the  thrust  alone.  A  blow  down  from 
00  =  50  to  00  =  20  is  replotted  in  Fig.  2  where  the 
curve  for  0o  =  20  is  normalized  by  the  same  pressure 
as  for  00  =  50 ,  making  all  curves  proportional  to  the 
thrust.  This  is  accomplished  by  multiplying  the  thrust 
coefficient  of  Fig  1 .  by  0.4  for  0o  =  20 .  The  curve  la¬ 
beled  “b.d.”  gives  the  average  thrust  coefficient  for  the 
blow  down  and  is  seen  to  be  an  average  of  the  curves 
for  00  =  50  and  0o  =  20 ,  as  expected.  The  variable 
nozzle  case  shown  in  Fig.  2  will  be  discussed  later. 

Average  thrust  coefficients  are  plotted  as  functions 

of  00  £  for  blow  downs  to  r  =  0.75 , 0.5,  and  0.25 

in  Figs.  3,  4,  and  5,  respectively.  Dimensionless  blow 
down  times  are  plotted  as  a  function  of  r  for  two  values 
of  yin  Fig.  6.  As  can  be  seen  from  eq.  (11)  and  Fig.  6, 
blow  down  times  are  independent  of  0o ,  and  depend 
only  weakly  on  y.  Because  t  is  fixed  and  independent  of 
00  for  a  fixed  r.  the  curves  in  Figs.  3,  4,  and  5  will  also 

be  proportional  to  the  dimensionless  impulse  £,  as  can 
be  shown  from  eq.  (13).  However,  the  constant  of  pro¬ 
portionality  will  be  different  for  each  r,  increasing  as  r 
decreases,  because  T  increases  as  r  decreases.  The  con¬ 
stant  of  proportionality  is  given  in  the  figures.  The 
overall  trend  in  these  figures  is  that  the  average  thrust 
coefficient  decreases  as  r  decreases.  This  is  due  to  t 
increasing  more  rapidly  than  i  as  r  decreases,  because  i 
also  increases  as  r  decreases. 

The  expansion  ratio  which  produces  the  maximum 
impulse  can  in  principle  be  found  by  taking  the  first 
derivative  of  eq.  (10)  with  respect  to  and  setting  it 
equal  to  zero.  The  resulting  expression  cannot  be  solved 
analytically,  however,  so  the  optimum  expansion  ratio 
is  computed  here  by  finding  the  maximum  impulse  us¬ 
ing  a  numerical  search.  The  optimum  expansion  ratio  is 
plotted  as  a  function  of  0o  and  r  in  Fig.  7,  and  the  di¬ 
mensionless  impulse  produced  at  the  optimum  expan¬ 
sion  ratio  is  plotted  in  Fig.  8.  Because  r  is  variable,  ris 
also,  so  the  dimensionless  impulse  and  the  average 
thrust  coefficient  curves  are  no  longer  directly  propor¬ 


tional  as  they  were  in  earlier  figures,  where  r  was  fixed. 
The  optimum  thrust  coefficient  is  plotted  in  Fig.  9. 

3.3  Variable  Nozzles 

Because  the  expansion  ratio  is  always  optimized 
for  variable  nozzles,  the  blow  down  of  variable  nozzles 
proceeds  along  the  curve  of  maximum  thrust  in  Fig.  1. 
The  dimensionless  impulse  produced  can  be  computed 
as  a  function  of  0o  and  r  by  integrating  eq.  (14),  and 
compared  with  the  performance  of  fixed  nozzles  oper¬ 
ating  at  the  optimum  expansion  ratio.  The  results  are 
not  easily  viewed  in  the  form  of  Fig.  8,  however,  so 
they  are  plotted  instead  here  in  Fig.  10  as  the  impulse 
penalty  in  using  optimized  fixed  nozzles  compared  to 
the  variable  nozzle  case,  based  on  the  percentage  of  the 
variable  nozzle  impulse.  The  penalty  is  shown  not  to 
exceed  about  3%  for  all  cases  calculated.  The  reason 
the  penalty  is  not  more  severe  can  be  traced  to  the  flat¬ 
ness  of  the  curves  in  Fig.  1  near  the  points  of  maximum 
thrust.  This  is  illustrated  more  clearly  in  Fig.  2,  where 
the  difference  between  the  variable  and  optimized  fixed 
nozzle  blow  down  paths  is  shown.  As  can  be  seen,  the 
maximum  instantaneous  thrust  coefficient  is  not  that 
much  higher  for  the  variable  nozzle  case  than  for  opti¬ 
mized  fixed  nozzles,  leading  to  calculated  average 
thrust  coefficients  which  are  very  close. 

3.4  Comparison  with  Constant  Pressure  Devices 

By  normalizing  with  the  same  constants  as  used  for 
the  CV  device,  the  correction  factors  Q  and  0  as  de¬ 
fined  in  eq.  (23)  and  used  in  eqs.  (20)-(22)  allow  the 
normalized  quantities  for  CP  devices  to  be  directly 
compared  with  CV  devices  at  any  arbitrary  CP  chamber 
pressure  and  speed  of  sound  (temperature).  However, 
only  certain  cases  will  be  of  interest  here.  It  is  pre¬ 
sumed  that  the  two  cycles  are  to  be  compared  for  the 
same  propellant  combination  and  at  the  optimum  mix¬ 
ture  ratio  for  each  cycle,  the  latter  of  which  as  discussed 
above  should  be  approximately  the  same  in  both  cases. 
If  Tf  is  the  initial  fill  temperature  before  combustion, 

and  q  is  the  heat  of  reaction  per  unit  mass  for  a  given 
propellant  combination  and  mixture  ratio,  then  the  ini¬ 
tial  CP  temperature  after  combustion  will  be  approxi¬ 
mately  T^p=Tf  +  q!  Cp  ,  while  the  initial  CV  tempera¬ 
ture  after  combustion  will  be  approximately 
Tq  -  Tf  +'Yq ICp  .  Thus  £2  will  vary  only  over  a  very 

narrow  range  near  unity,  namely  from  Q  -  \  for  q  =  0 
to  Q  for  qlCp  »Tf  .  For  y  =1.2,  this  corre¬ 
sponds  to  1  <  <  1 .095  . 

The  next  question  is  which  of  the  various  pressures 
encountered  over  the  CV  blow  down,  or  not  encoun¬ 
tered,  should  be  used  for  a  CP  comparison.  Only  two 
possibilities  will  be  considered  here.  The  first  is  the 
case  where  the  pressure  of  the  CP  device  is  the  same  as 
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the  fill  pressure  of  the  CV  device,  namely  0  =  Pq  I . 

This  would  allow  a  comparison  of  two  devices  having 
equivalent  feed  systems.  Assuming  the  purge  pressure 
is  the  same  as  the  feed  pressure,  limit  cycle  operation 

gives  r  =  (Pj- / .  Given  further  that  Pq  =  PqRTq 

and  that  Pg  =  Pf  for  a  CV  device,  where  Py  is  the 

density  after  propellant  fill  but  before  combustion,  then 
PqIPj  =10  would  be  representative  of  propellants  in 

an  initial  gaseous  state,  while  Pq! Pf  =1000  would 

approach  propellants  in  an  initial  liquid  state.  A  CP 
device  operating  at  the  same  pressure  as  the  fill  pres¬ 
sure  of  the  CV  device  will  be  referred  to  here  as  a  “CP/’ 
device  for  brevity. 

The  second  possibility  will  be  to  compare  the  per¬ 
formance  of  combustion  chambers  experiencing  the 
same  material  stresses.  A  very  rough  approximation  of 
this  would  be  a  CP  device  operating  at  the  same  pres¬ 
sure  as  the  peak  pressure  of  a  CV  device,  namely 
0=1.  A  CP  device  operating  at  the  same  pressure  as 
the  peak  pressure  of  the  CV  device  will  be  referred  to 
here  as  a  “CPp”  device,  again  for  brevity. 

The  ratio  of  the  dimensionless  impulse  produced 
by  a  CPy  device  to  the  dimensionless  impulse  of  a  CV 
device  is  presented  in  Fig.  1 1 ,  and  the  ratio  of  the  di¬ 
mensionless  impulse  produced  by  the  CPp  device  to  the 
dimensionless  impulse  of  a  CV  device  is  presented  in 
Fig.  12.  Because  the  same  quantities  are  used  to  nondi- 
mensionalize  the  impulse  in  all  cases,  these  plots  also 
give  the  ratio  of  the  total  impulses  and  the  ratio  of  the 
specific  impulses.  The  expansion  ratios  used  for  the  CP 
and  CV  devices  are  those  which  produce  the  maximum 
possible  impulse  in  each  case,  namely  a  fixed  expan¬ 
sion  ratio  which  matches  the  fixed  pressure  ratio  for  the 
CP  device,  and  a  variable  expansion  ratio  to  continu¬ 
ously  match  the  variable  pressure  ratio  for  the  CV  de¬ 
vice.  The  CV  device  is  operated  in  limit  cycle  mode 
where  the  blow  down  is  performed  to  the  maximum 
value  of  r  consistent  with  limit  cycle  operation, 

r  =  {Pf !  Pof'^  .  It  is  assumed  that  qlCp»Tf  in  all 
cases,  giving  Q  =  ■ 

Curves  are  drawn  in  Figs.  1 1  and  12  for  different 
ratios  of  the  peak  pressure  to  the  fill  pressure,  P^l  Pj  . 

Increases  in  this  ratio  indicate  either  more  energetic 
chemistry  (greater  temperature/pressure  rise),  or  the 
presence  of  condensed  phases  in  the  initial  fill.  In¬ 
creases  over  two  orders  of  magnitude  as  shown  in  the 
figures  would  be  largely  the  result  of  the  presence  of 
condensed  phases.  Increases  in  Pq  /  Pf  always  increase 

the  difference  between  the  CP  and  CV  performance. 
The  right  hand  extent  of  all  the  curves  is  dictated  by  the 
shock  limit,  beyond  which  a  shock  enters  the  nozzle 


and  the  present  formulation  no  longer  applies.  The  left 
hand  extent  of  the  curves  as  drawn  in  the  figures  is  lim¬ 
ited  solely  for  reasons  of  graphical  clarity;  the  curves  in 
both  figures  actually  all  extend  to  l/0o  =  0  •  The  limit¬ 
ing  values  of  the  curves  at  \I^q=0  are  given  in  the 
figures.  These  limiting  values  are  exactly  the  same  for 
both  the  CPy  and  CPp  cases.  The  reason  for  this  is  that 
eq.  (20)  becomes  independent  of  the  chamber  pressure 
when  1/^0  =  0  ■ 

Away  from  1/(I)q  =  0 ,  the  impulse  of  a  CP/  device 
is  shown  in  Fig.  11  to  be  inferior  to  that  of  the  CV  de¬ 
vice,  and  the  impulse  of  a  CPp  device  is  shown  in  Fig. 
12  to  be  superior  to  that  of  the  CV  device.  The  relative 
differences  diminish  as  1/ 0q  approaches  zero.  At  large 
enough  l/^g.  relative  differences  are  significantly 
larger  than  the  3%  maximum  performance  penalty 
shown  in  Fig.  10  attributed  to  using  an  optimized  fixed 
nozzle  instead  of  a  variable  nozzle.  This  implies  that,  at 
large  enough  l/0o  -  the  performance  penalty  associated 
with  the  practical  simplicity  of  using  an  optimized  fixed 
nozzle  might  be  tolerable  compared  with  the  complica¬ 
tions  of  developing  a  variable  nozzle.  To  assess 
whether  a  practically  attainable  device  could  fall  within 
a  regime  where  such  losses  might  be  tolerable,  a  refer¬ 
ence  case  consisting  of  the  approximate  operating  point 
of  a  stoichiometric  kerosine/air  (jet  fuel/air)  system 
operating  at  1  atm  ambient  pressure  and  a  3  atm  fill 
pressure  is  plotted  as  a  black  dot  in  Fig.  11.  The  oper¬ 
ating  point  is  shown  to  be  well  within  a  regime  where 
such  performance  losses  might  be  tolerable.  Whether 
such  losses  would  in  fact  be  tolerable  in  a  practical 
system  would  of  course  depend  on  all  the  other  system 
losses. 

The  case  l/0o  =  0  important  limit  in  that  it 
represents  operation  where  the  ambient  pressure  is  a 
vacuum,  i.e.,  in  space.  As  l/0o  approaches  zero,  the 
difference  between  a  CP/  device  and  a  CPp  device  van¬ 
ishes,  and  the  impulse  of  the  CP  device  can  become 
either  slightly  higher  or  slightly  lower  than  the  impulse 
of  the  CV  device,  depending  on  the  magnitude  of 
PqI  Pf .  The  reason  for  this  may  be  found  in  deriving 

the  following  analytical  expression  for  i/i^p,  which  is 
possible  when  l/0o  =0.  For  limit  cycle  operation,  in¬ 
tegrate  eq.  (14)  from  r  =  {PjlPQf'''  to  r  =  l,  which 

now  can  be  done  analytically,  and  add  the  impulse  from 
the  constant  pressure  part  of  the  cycle  by  using  the 
functions  from  eq.  (18)  in  eq.  (10).  Note  that  =0 
(£  =  oo )  when  1/00  =  0  •  Divide  the  result  into  eq.  (20). 
The  result  is 
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The  maximum  value  of  the  expression  occurs  as 
which  for  Q  and  y  =1.2  is 

ili^p  =1.004.  As  PqIPj  becomes  smaller,  the  magni¬ 
tude  of  the  impulse  ratio  also  becomes  smaller,  until 
eventually  it  is  reduced  to  a  value  less  than  unity.  For 
Po^Pf  =10,  this  value  becomes  i/i^p  =0.992.  Limit 

values  of  ili^p  at  l/0o  shown  in  Figs  11  and  12 
were  calculated  using  eq.  (24). 

4.0  Conclusions 

Analytical  solutions  of  the  constant  volume  limit  of 
pulsed  propulsion  for  fixed  expansion  ratios  have  been 
explored  and  compared  with  the  case  where  a  variable 
nozzle  is  used  to  match  the  pressure  ratio  at  all  times 
during  the  blow  down.  The  solutions  apply  for  all  su¬ 
personic  flows  where  compression  waves  remain  exte¬ 
rior  to  the  nozzle.  It  was  found  that: 

1.  CV  devices  should  optimize  at  approximately  the 
same  mixture  ratio  as  CP  devices,  namely  the  mix¬ 
ture  ratio  which  maximizes  the  initial  temperature 
and  minimizes  the  molecular  weight. 

2.  The  thrust  of  CV  devices  is  maximized  in  the  same 
way  as  for  CP  devices,  namely  by  maximizing  the 
initial  pressure  and  the  throat  area. 

3.  The  blow  down  time  depends  on  r  =  p/po,  only 
weakly  on  y,  and  is  independent  of  the  initial  pres¬ 
sure  ratio  <I>Q  =  Pq/  P„  . 

4.  In  general  an  optimum  expansion  ratio  exists  which 
maximizes  the  impulse  produced  by  the  blow  down 
of  a  CV  device  at  a  fixed  expansion  ratio. 

5.  Using  an  optimized  fixed  expansion  ratio  nozzle 
results  in  an  impulse  penalty  that  does  not  exceed 
3%  (for  the  cases  considered)  of  the  impulse  that 
would  be  produced  using  the  more  complicated 
variable  expansion  ratio  nozzle. 

6.  A  comparison  between  a  CP  device  operating  at  its 
optimum  expansion  ratio  and  a  CV  device  in  limit 
cycle  operation  with  a  variable  expansion  ratio  leads 
to  the  following  additional  conclusions: 

a.  Except  near  1/0q  =  F^/Pg  =0  ,  the  impulse 
produced  by  a  CV  device  is  superior  to  that  of  a 
CP  device  operating  at  the  same  pressure  as  the 
fill  pressure  of  the  CV  device,  but  inferior  to  the 
impulse  produced  by  a  CP  device  operating  at 
the  same  pressure  as  the  peak  pressure  of  the 
CV  device. 

b.  The  magnitude  of  the  difference  between  a  CV 
device  and  the  two  CP  devices  considered  in¬ 
creases  as  1/00  increases. 


c.  At  large  enough  1/0O'  relative  difference 
between  a  CV  device  and  the  two  CP  devices 
considered  is  significantly  larger  than  the  im¬ 
pulse  penalty  associated  with  using  an  opti¬ 
mized  fixed  expansion  ratio  instead  of  a  vari¬ 
able  nozzle  on  the  CV  device. 

d.  The  regime  where  the  relative  difference  be¬ 
tween  a  CV  device  and  a  CP  device  is  signifi¬ 
cantly  larger  than  the  impulse  penalty  associ¬ 
ated  with  using  an  optimized  fixed  expansion 
ratio  on  the  CV  device  appears  to  be  practically 
achievable  for  at  least  the  one  reference  case 
considered,  which  involved  kerosine/air  opera¬ 
tion  at  1  atm  ambient  pressure. 

e.  Near  1/0q  =  P^/Pq  =  0  ,  which  is  an  important 
limit  in  that  it  represents  operation  at  the  near 
zero  ambient  pressure  of  space,  the  difference 
between  the  two  CP  devices  considered  van¬ 
ishes,  and  the  impulse  of  the  CP  device  can  be¬ 
come  either  slightly  higher  or  slightly  lower 
than  the  impulse  of  the  CV  device,  depending 
on  the  magnitude  of  Pq/  P/  ■ 
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cerning  this  work  are  gratefully  acknowledged. 

References 

1.  Kailasanath,  K.,  Patnaik,  G.,  and  Li,  C., 

“Computational  Studies  of  Pulse  Detonation 
Engines:  A  Status  Report,”  35'*' 

AIAA/ASME/SAE/ASEE  Joint  Propulsion 
Conference  and  Exhibit,  paper  AIAA  99-2634,  20- 
24  June  1999,  Los  Angeles,  CA. 

2.  Leipman,  H.W.,  and  Roshko,  A.,  Elements  of 
Gasdynamics,  John  Wylie  &  Sons,  1957. 

3.  Sutton,  G.P.,  Rocket  Propulsion  Elements,  6'*'  ed., 
John  Wiley  &  Sons,  1992. 


7 


Figure  1.  Steady  state  thust  coefficients 


Figure  2.  Blow  down  from  =  50  to  (/>o  =  20  ,  where 
“b.d.”  curve  is  for  the  blow  down.  Arrows  indicate  paths 
taken  for  fixed  and  variable  nozzles. 
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Figure  5.  Blow  down  to  r  =  0.25 


Figure  6.  Blow  down  times 


Figure  7.  Optimum  expansion  ratios 


Figure  8.  Dimensionless  impulse  at  the  optimimum  expansion  ratio 
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Figure  9.  Thrust  coefficient  at  the  optimum  expansion  ratio 


Figure  11.  Comparison  of  CP/  and  CV  impulse.  The  reference 
case  is  stoichiometric  kerosine/air  at  1  atm  ambient  pressure  and 
3  atm  fill  pressure. 
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Figure  10.  Impulse  penalty  of  an  optimized  fixed  nozzle 
Compared  with  a  variable  nozzle 


Figure  12.  Comparison  of  CP^  and  CV  impulse. 
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